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We show that partial dynamical symmetries (PDS) can occur at critical-points of quantum phase 
transitions, in which case, underlying competing symmetries are conserved exactly by a subset of 
states, and mix strongly in other states. Several types of PDS are demonstrated with the example of 
critical-point Hamiltonians for first- and second-order transitions in the framework of the interacting 
boson model, whose dynamical symmetries correspond to different shape-phases in nuclei. 
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Symmetries play a profound role in determining the 
critical behaviour of dynamical systems. Their signifi- 
cance was recognized in Landau's classic theory of ther- 
mal phase transitions [l[ and in the renormalization 
group of critical phenomena Q. An equally- important 
role is played by symmetries in quantum phase transi- 
tions (QPT) or ground-state energy phase transitions 
which occur at zero temperature as a function of a cou- 
pling constant. Such structural changes are currently 
of great interest in different branches of physics. QPT 
occur as a result of a competition between terms in the 
Hamiltonian with different symmetries which lead to con- 
siderable mixing in the eigenfunctions, especially at the 
critical-point where the structure changes most rapidly. 
In the present work we address the question whether 
there are any symmetries (or traces of) still present at the 
critical-points of QPT. We show that particular symme- 
try constructions, called partial dynamical symmetries, 
can survive at the critical-point in-spitc of the strong 
mixing. The feasibility of such persisting symmetries 
gains support from the recently proposed |4| and empir- 
ically confirmed [5] analytic descriptions of critical-point 
nuclei, and the emergence of "quasi-dynamical symme- 
tries" @ in the vicinity of such critical-points. 

A convenient framework to study symmetry-aspects of 
QPT are models where the Hamiltonian is expanded in 
elements of a Lie algebra (Go), called the spectrum gen- 
erating algebra. A dynamical symmetry occurs if the 
Hamiltonian can be written in terms of the Casimir op- 
erators of a chain of nested algebras of Go, 



Go D G\ D ■ ■ ■ D G n 



(1) 



terminating with an invariant algbera G n . The following 
properties are then observed, (i) All states are solvable 
and analytic expressions are available for energies and 
other observables. (ii) All states are classified by quan- 
tum numbers, |ao, oli , ■ ■ • , at n ), which are the labels of the 
irreducible representations (irreps) of the algebras in the 
chain, (iii) The structure of wave functions is completely 
dictated by symmetry and is independent of the Hamil- 
tonian's parameters. Partial dynamical symmetry (PDS) 
corresponds to a particular symmetry-breaking for which 
some (but not all) of the above mentioned virtues of a 



dynamical symmetry are retained. PDS of type I corre- 
sponds to a situation where some of the states have all 
the dynamical symmetry. In this case the properties (i)- 
(iii) are fulfilled exactly, but by only a subset of states. 
PDS of type II corresponds to a situation for which all the 
states preserve part of the dynamical symmetry. In this 
case there are no analytic solutions, yet selected quan- 
tum numbers (of the conserved symmetries) are retained. 
This can occur, for example, when, the Hamiltonian pre- 
serves only selected symmetries Gi C G n in the chain (JTJ) , 
and only their irreps are unmixed. PDS of type III has 
a hybrid character, for which some of the states preserve 
part of the dynamical symmetry. PDS of various types 
have been shown to be relevant to nuclear and molecular 
spectroscopy 0, 0, S E3, EE E3 and to mixed systems 
with coexisting regularity and chaos All examples 

of PDS encountered so far involved stable limits of struc- 
ture. In the present work we show the relevance of the 
PDS notion to the more complicated case of a phase tran- 
sition. 

As a concrete example, we consider the interacting bo- 
son model (IBM) widely used in the description of 
quadrupole collective states in nuclei, in terms of a sys- 
tem of N monopole (s) and quadrupole (al) bosons, rep- 
resenting valence nucleon pairs. The spectrum generat- 
ing algebra is Go = U (6) and the invariant algebra is 
G n = 0(3). The three dynamical symmetry limits of the 
model and corresponding bases are 



17(6) D 17(5) D 0(5) D 0(3) 
17(6) D SU{3) D 0(3) 
17(6) D 0(6) D 0(5) D 0(3) 



\N,n d ,r,v,L) (2a) 
\N,(X,n),K,L) (2b) 
\N,a,T,P,L) . (2c) 



The quantum numbers N, rid, (A, fj,),a,T and L label the 
relevant irreps of C/(6), 17(5), SI7(3), 0(6), 0(5) and 0(3) 
respectively, v and K are multiplicity-labels needed for 
complete classification of selected states in the reductions 
0(5) D 0(3) and 577(3) D 0(3) respectively. The ana- 
lytic solutions of these dynamical symmetries resemble a 
spherical vibrator, axially-deformed rotor and deformed 
7-soft rotor for the C7(5), SU(3) and 0(6) chains respec- 
tively. This identification is consistent with the geometric 
visualization of the model in terms of a potential surface 
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FIG. 1: (color online). Spectrum of H(/3 = y/2), Eq. ©, 
with h 2 = 0.1 and N = 10. L = 2 ,3i are solvable 17(5) 
states of Eq. {9). L(A" = 0i) and L(if = 2i) are, respectively, 
solvable 517(3) states of Eq. (10a) and Eq. (10b) with k = 1. 
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FIG. 2: U(5) (rid) and 5(7(3) [(A, fi)] decomposition of se- 
lected spherical and deformed states in Fig. 1. 



defined by the expectation value of the Hamiltonian in 
the coherent (intrinsic) state [H, [l6| 



|/3, 7 ; AO = (N\)-V 2 (bt) N \0) 



(3) 



where b\ = {l+f3 2 )- 1 / 2 [f3cos 1 dl+(3sm 1 {d\+d ] _ 2 )/^+ 
s T ]. For the general IBM Hamiltonian with one- and two- 
body interactions, the potential surface reads 



E(J3,i)=Eo + N(N-l) 



6/3 3 cos 3 7 + c/3 4 



(1 + /3 2 ) 2 



•(4) 



The coefficients Eo,a,b,c involve particular linear com- 
binations of the Hamiltonian's parameters 0. The 
quadrupole shape parameters 7) at the global min- 
imum of E((3, 7) define the equilibrium shape for a given 
Hamiltonian. The shape can be spherical ((3 = 0) or de- 
formed (/3 > 0) with 7 = (prolate), 7 = n/3 (oblate), 
or 7-independent (6 = 0). 

Phase transitions can be studied by IBM Hamiltonians 
of the form, H± + gH 2 , involving terms from different dy- 
namical symmetry chains flij ]. The nature of the phase 
transition is governed by the topology of the correspond- 
ing surface ([4|) , which serves as a Landau's potential with 
the equilibrium deformations as order parameters. The 
conditions on the surface at the critical-points of first- 
and second-order transitions are 



I s * order 
2 nd order 



b 2 = 4ac, a > 0, 6^0 
o = 0, 6 = 0, c>0. 



(5a) 
(5b) 



The first-order critical-surface (5a) has degenerate spher- 
ical and deformed minima at (3 = and ((3 — 2<z/|6|,7o), 
with 70 = 0(7r/3) for b > 0(6 < 0). The second-order 
critical-surface (5b) is independent of 7 and behaves as 
(3 A for small [3. The conditions in Eq. (5) fix the critical 
value of the control parameter (g — g c ) which, in turn, 



determines the critical-point Hamiltonian. IBM Hamilto- 
nians of this type have been used extensively for studying 
shape-phase transitions in nuclei [il ITj. 



19, 20, 21 . We 



now show that a large class of such Hamiltonians exhibit 
PDS. 

The spherical to deformed 7-soft shape-phase transi- 
tion is modeled in the IBM by the Hamiltonian 



H = 



4(JV 



■1)A, 



(^) 2 ][H.c] 



(6) 



where H.c. stands for Hermitian conjugate and the dot 
implies a scalar product. The n^-term is the d-boson 
number operator (eigenvalues rid), which is the linear 
Casimir operator of U(5). The A-teim is related to the 
Casimir operator of 0(6) (HI- For the indicated ratio of 
coefficients, the above H satisfies condition (5b), hence 
qualifies as a second-order critical Hamiltonian. The first 
(second) term in Eq. ((6]) has 0(6) [f/(5)] selection rules 
Act = 0,±2 (An d = 0,±2), and both terms are 0(5)- 
scalars. Consequently, the eigenstates of H have good 
0(5) symmetry (r), but are mixed strongly with respect 
to both (7(5) and 0(6) 0. Since both (7(5) and 0(6) 
are broken while 0(5) D 0(3) are preserved, by defini- 
tion, the critical Hamiltonian has an 0(5) PDS of type II. 
In fact, since 0(5) is a good symmetry common to both 
chains (2a) and (2c), the 0(5) PDS is valid throughout 
the J7(5)-0(6) transition region. 

A recent study of QPT within the IBM has shown 
that, apart from rotational terms which do not affect the 
potential surface of Eq. the critical Hamiltonian for a 
spherical to prolate-deformed shape-phase transition can 
be transcribed in the form [211 



H(/3 ) = h 2 P^{/3 ) ■ P2((3 ) 



(7) 



where £*,(/%) = ftatrft + ^7/2 (dU^\ P 2 ,(f3 ) = 
(— iy_P 2 ,-/i(/3o) and h 2 , A) > 0- The corresponding sur- 
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face in Eq. ([4]) has coefficients a = h 2 0Q, b = 2h 2 (3 , c = 
h 2 , which satisfy condition (5a). This qualifies H((3 ) as 
a first-order critical Hamiltonian whose potential accom- 
modates two degenerate minima at f3 — and (/3, 7) = 
(/3a, 0). H(/3o) is constructed to have the equilibrium 
intrinsic state, |/3 = Poil — 0;N), Eq. $5§, as a zero- 
energy eigenstate. Rotational-invariance ensures that 
states, \(3q;N,L), of good 0(3) symmetry L projected 
from this intrinsic state, remain zero-energy eigenstates. 
H((3o) then has a solvable deformed ground band, 

W,N,L) E = Q (L = 0,2,4,..., 27V) . (8) 

It has also the following solvable spherical eigenstates 

\N,n d = t = L = 0) E = (9a) 
\N,n d = t = L = 3) E = 3h 2 \Pl(N -3) + S\ . (9b) 

As shown in Figs. 1-4, the remaining states in the spec- 
trum of H(Po) are either predominantly spherical (with 
characteristic dominance of single rid components) or de- 
formed states (with a broad rid distribution) arranged in 
several excited bands. 

The critical Hamiltonian of Eq. |(7J) with (3q = \[2 is a 
special case of a Hamiltonian shown in 0] to have SU(3) 
PDS. This comes about because the sequence of states 
l fc > « [ P 2,2(V2)} k \f3 = y/2, 7 = 0; N - 2k) are eigenstates 
of H{flo = y/2). These are lowest-weight states in the 
SU(3) irreps (A, fi) = {2N - 4k, 2k) with 2k < N. In the 
nuclear physics terminology they are referred to as in- 
trinsic states representing deformed ground (k — 0) and 
7 fe bands, with angular momentum projection (K = 2k) 
along the symmetry axis. Since H([3q — y/2) is an 0(3)- 
scalar, the states of good L projected from \k) remain 
eigenstates with quantum numbers, \N, (A, n)K, L), of 



the SU(3) chain (2b), and form solvable bands, 

|JV, (2N, 0)K = 0, L) E = [L = 0,2,4, ... ,2N) (10a) 
\N, (2N - 4k, 2k)K = 2k, L) E = 3h 2 [2N + 1 - 2k]k 
L = K,K + 1,K + 2,...,(2N -2k) k > . (10b) 

In addition, H((3 — v2) has the spherical states of 
Eq. ©, with good U(5) symmetry, as eigenstates. The 
remaining levels of H(/3q = V2), shown in Fig. 1, are 
calculated numerically. Their wave functions are spread 
over many U(5) and SU(3) irreps, as is evident from 
Fig. 2. This situation, where some states are solvable 
with good t/(5) symmetry, some are solvable with good 
SU(3) symmetry and all other states are mixed with re- 
spect to both U(5) and SU(3), defines a C/(5) PDS of 
type I coexisting with a SU(3) PDS of type I. 

The Hamiltonian of Eq. |[7J) with (3n = 1 is a spe- 
cial case of a Hamiltonian shown in [9j to have 0(6) 
PDS. This comes about because the intrinsic state of 
Eq. ([3|) with (0o = 1,7 = 0) is a zero-energy eigenstate 
of H(P = 1) with good 0(6) symmetry (a = N). The 
0(3)-invariance of the Hamiltonian ensures that states of 
good L projected from \/3q = 1,7 = 0; N) form a solvable 
ground band with good 0(6) character, 

\N,a = N,L) E = (L = 0, 2, 4, . . . , 2N) . (11) 

In addition, if (/3q = 1) has the spherical states of Eq. |9]), 
with good U (5) symmetry, as eigenstates. The remaining 
eigenstates in Fig. 3 are mixed with respect to both U(5) 
and 0(6), as is evident from their decomposition shown 
in Fig. 4. Apart from the solvable U(b) states of Eq. ((9]), 
all eigenstates of H(/3q = 1) are mixed with respect to 
0(5) (including the solvable 0(6) states of Eq. (fTTjl . as 
shown in Fig. 5). It follows that the Hamiltonian has a 
subset of states with good U(5) symmetry and a subset 
of states with good 0(6) but broken 0(5) symmetry, and 
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FIG. 3: (color online). Spectrum of H(/3 = 1), Eq. ©, with 
hz = 0.1 and N = 10. L = 02,3i are solvable U(5) states of 
Eq. ©. L(K — 0i) are solvable states of Eq. |TT} with good 
0(6) but broken 0(5) symmetry. 
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FIG. 4: U(5) (rid) and 0(6) (a) decomposition of selected 
spherical and deformed states in Fig. 3. 
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FIG. 5: 0(5) (r) decomposition of the L = 0, 2 states, 
Eq. members of the ground band (A" = 0i) shown in 

Fig. 3. Both states have 0(6) symmetry a = N. 



all other states are mixed with respect to both Z7(5) and 
0(6). These are precisely the required features of £7(5) 
PDS of type I coexisting with 0(6) PDS of type III. 

For arbitrary values of 0o, the spherical states with 
good £7(5) symmetry, Eq. (|5J), are still eigenstates of 
H(p ), Eq. 0, hence £7(5) PDS of type I is still valid. In 
general, the deformed states of Eq. |8|), are not associated 
with any IBM dynamical symmetry but, nevertheless, are 
still solvable. This situation may be referred to as partial 
solvability. Since the wave functions of the solvable states 
are known, it is possible to obtain closed form expressions 
for related observables. For example, for the electromag- 
netic E2 operator, T{E2) = Ss + d + X (rf d) {2 \ the 
necessary matrix elements for transitions involving states 
inEqs. ©-© are Ti = (/?; N, L\ \ T(E2) \\/3;N,L+2) and 
T 2 = (/?; N,L = 2\\T(E2) \\N,n d = r = L = 0), 



fla 1 rffl 1 (/3) + a a r£r*'Q9) 

C9)]Va 



Ti 
T 2 



Cr 



(L+2) 



(12) 



where C L = ^2L + 5(L + 2,0;2,0|L,0) is proportional 
to a Clebsch Gordan coefficient, a\ = 1 — PxL/ (2-£ + 3), 
a 2 = l-/3x(L+3)/(2L+3) with * = v^X, an d r^(/3) 
is a normalization factor given in [2 11 ] - 

As discussed, the spectrum of H(f3 ), Eq. exhibits 
coexistence of spherical and deformed states, signaling a 
first-order transition. In particular, the spherical L = 
state, Eq. (9a), is exactly degenerate with the ground 
band, Eq. (J8j) , and for /3 = \[2 also the spherical L = 3 
state, Eq. (9b), is degenerate with the SU(3) 7-band, 
Eq. (10b) with k = 1. Adding to the Hamiltonian the 
Casimir operator of 0(3), contributes an exact L(L + 1) 
splitting with no effect on wave functions. The remaining 
degeneracy of states with the same L, can be lifted by 
adding a small one-body term rid- With that, the spher- 
ical U(5) states of Eq. §§§ remain solvable eigenstates. 
However, the fid term destroys the exact solvability and 



partial-symmetry of the deformed states, Eq. |8]). The 
corresponding leading-order shifts can be estimated from 
(/?; N, L\n d \f3; N,L)=N- T^P) /T^ ((3) . 

In summary, we have shown the relevance of the PDS 
notion to critical-points of QPT, with phases character- 
ized by Lie-algebraic symmetries. In the example consid- 
ered, second-order critical Hamiltonians mix incompati- 
ble symmetries but preserve a common lower symmetry, 
resulting in a single PDS with selected quantum num- 
bers conserved. First-order critical Hamiltonians exhibit 
distinct subsets of solvable states with good symmetries, 
giving rise to a coexistence of different PDS. The ingre- 
dients of an algebraic description of QPT is a spectrum 
generating algebra and an associated geometric space, 
formulated in terms of coherent (intrinsic) states. The 
same ingredients are used in the construction of Hamilto- 
nians with PDS. These, in accord with the present work, 
can be used as tools to explore the role of possibly partial 
symmetries in governing the critical behaviour of diverse 
dynamical systems undergoing QPT. This work was sup- 
ported by the Israel Science Foundation. 
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